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Abstract 

Background: When diseased intestine (i.e., from colon cancer, diverticulitis) requires resection, its reconnection 
(termed anastomosis) can be complicated by non-healing of the newly joined intestine resulting in spillage of 
intestinal contents into the abdominal cavity (termed anastomotic leakage). While it is suspected that the intestinal 
microbiota have the capacity to both accelerate and complicate anastomotic healing, the associated genotypes and 
functions have not been characterized. 

Results: Using 16S rRNA amplicon sequencing of samples collected on the day of surgery (postoperative day 0 
(PODO)) and the 6th day following surgery (postoperative day 0 (P0D6)), we analyzed the changes in luminal versus 
tissue-associated microbiota at anastomotic sites created in the colon of rats. Results indicated that anastomotic injury 
induced significant changes in the anastomotic tissue-associated microbiota with minimal differences in the luminal 
microbiota. The most striking difference was a 500-fold and 200-fold increase in the relative abundance of Enterococcus 
and Escherichia/Shigella, respectively. Functional profiling predicted the predominance of bacterial virulence-associated 
pathways in post-anastomotic tissues, including production of hemolysin, cytolethal toxins, fimbriae, invasins, cytotoxic 
necrotizing factors, and coccolysin. 

Conclusion: Taken together, our results suggest that compositional and functional changes accompany anastomotic 
tissues and may potentially accelerate or complicate anastomotic healing. 

Keywords: Colon anastomosis, 16S rRNA, PiCRUST, Bacterial composition, Predicted function, Anastomotic tissues, 
Luminal content 



Background 

Removal (resection) of a segment of the gastrointestinal 
tract is a common procedure performed by general sur- 
geons for a variety of illnesses including cancer, diver- 
ticulitis, or intestinal obstruction. When surgeons resect 
a diseased segment of the gastrointestinal tract, intes- 
tinal continuity is reestablished by carefully suturing or 
stapling the remaining ends together to create a viable 
and well-sealed connection referred to as an 'anasto- 
mosis'. An anastomosis therefore is a significant tissue 
injury given that the surgeon must first divide all of the 
blood vessels feeding the proposed resection site; all 
layers of the intestinal wall are divided meaning that 
their entire thickness from mucosa to serosa is disrupted 
and then tissues are approximated with sutures and 
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staples that represent foreign bodies embedded into the 
area of healing tissues. As such, surgeons routinely ad- 
minister antibiotics directed against the intestinal micro- 
biota to decontaminate the intestine of its microflora in 
an effort to reduce the presumed effect of the microbiota 
on healing [1]. Despite excellent suture technique by 
highly skilled surgeons to approximate the tissues and 
fuse and seal the intestinal wall, without proper subse- 
quent healing of the anastomotic tissues, 'leak' can occur 
resulting in spillage of intraluminal contents into the 
peritoneal cavity that is associated with sepsis, permanent 
need for an ostomy, and death. Thus, despite excellent 
surgical technique, leaks still occur and its pathogenesis 
remains unknown. While surgeons have known for over 
60 years that the intestinal microbiota play a key and crit- 
ical role in the etiopathogensis of leakage [2], the species 
involved and the molecular details of their activity remain 
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unknown. Work from our lab has recently confirmed that 
both the species and phenotype of microbes present on 
anastomotic tissues is a critical determinant of anasto- 
motic leakage following intestinal resection and anasto- 
mosis [3], In this study, we identified that when high 
swarming phenotypes of Pseudomonas aeruginosa are 
present on anastomotic tissues, they are associated with 
leaks. Yet, the choice of antibiotics to provide prophylaxis 
to patients is based on culture-derived results of the lu- 
minal microbiota before and after antibiotic adminis- 
tration that are now several decades old. Because 
anastomotic tissues are generally inaccessible, it is ex- 
tremely difficult to track changes in microbial content 
over the course of healing. As a result, it is not known 
whether the community structure, membership, and/or 
function of the luminal microbiota are representative 
of the microbiota that associate with anastomotic tis- 
sues themselves. Furthermore, it is now recognized 
that culture-based studies of the luminal microbiota 
fail to classify the majority of the microbial diversity 
present. Finally, genetic sequencing has revealed the 
dependency of intestinal microbial community struc- 
ture, membership, and function on the regional and 
spatial context in which they establish niche specialization 
[4,5]. Therefore, there is a significant gap in knowledge as 
to whether healing anastomotic tissues harbor a distinct 
microbial community from the luminal microbial commu- 
nity. Elucidation of this issue may have a significant 
impact on the rationale for choosing antibiotics in 
preparation for colon surgery. 

Therefore, the aims of this study were to create anas- 
tomotic injury in rats similar to that which is performed 
in humans and to determine the changes and differences 
in the luminal versus tissue microbiota over a 6-day 
course of anastomotic healing. 



Results 

Anastomosis induces changes in the composition of the 
intestinal microbiota associated with intestinal tissue, but 
does not affect the microbiota associated with luminal 
contents 

Rats were subjected to removal (resection) of 1 cm of 
colon just above the colorectal junction followed by re- 
connection (anastomosis) as is typically performed in 
humans with colon cancer of the distal intestine. DNA 
was isolated from resected colon and colon luminal con- 
tent at the time of surgery, represented as postoperative 
day 0 (PODO). On postoperative day 6 (POD6), rats were 
sacrificed and the luminal and tissue associated micro- 
biota at the anastomotic segment were harvested for 
analysis. Results demonstrated that colon surgery did 
not affect the microbial community structure associated 
with the luminal contents (stool), but induced significant 
changes in the composition of the intestinal microbiota 
associated with the intestinal tissue (Figure 1A,B). 

The tissue microbiota at P0D6 is enriched by commensal 
opportunistic bacteria 

The taxonomic profiling of tissue-associated microbiota 
demonstrated that most of the changes in bacterial com- 
position were observed in the phyla Proteobacteria, Actino- 
bacteria, and Firmicutes (Figure 2). Among Proteobacteria, 
the relative abundance of Escherichia/Shigella {c_Gamma- 
proteobacteria;o_Enterobacteriales;i_Enterobacteriaceae) in- 
creased 200-fold. In Actinobacteria, the relative abundance 
of an uncultured bacterium (c_Coriobacteria;o_Corio- 
bacteriales;f_Coriobacteriaceae) increased approximately 
50-fold (Figure 3). Non-uniform changes among the genera 
Firmicutes were observed. The relative abundance of un- 
cultured bacteria representing Ruminococcaceae and Clos- 
tridia decreased 20-fold, while the relative abundance of 
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Figure 1 Principal components analysis of luminal- and tissue-associated microbiomes. Colon surgery affects the composition of the 
intestinal microbiota associated with intestinal tissue (B) but not the microbiota associated with luminal contents (A). n = 10 per group. 
PERMANOVA p values based on 10,000 random permutations of the dataset. Multiple comparisons were performed using the Bonferroni method. 
Error bars = ±2 standard error. 
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Figure 2 Taxonomic tree of anastomotic tissue microbiome changes. Shading indicates increases {green) and decreases (blue) in relative 
abundances of taxonomic groups on POD6 relative to PODO. Point sizes represent the higher point-in-time average abundance of the group 
relative to the entire microbial community. Groups comprising less than 1% of the total community on PODO and P0D6 are excluded. 



the genera Allobaculum and Coprococcus increased ap- 
proximately 10-fold, and Enterococcus increased 500-fold 
(Figure 3). Finally, a 20-fold decrease in the relative abun- 
dance of Prevotellaceae, demonstrated that Bacteroidetes 
were also influenced by anastomotic injury (Figure 3). 

Predicted functional profiling suggests the potential of 
virulence-associated pathways in post-anastomotic 
tissues 

The functional potential of the microbiota associated 
with the intraluminal contents (stool) and tissue before 
(PODO) and after anastomosis (POD6) was predicted 
using PiCRUST [6]. There were no statistically signifi- 
cant functional differences for the intraluminal contents 
between PODO and POD6, while the microbiota associ- 
ated with intestinal tissue revealed multiple significantly 
altered functional attributes after anastomosis (Table 1). 



The upregulated functions in tissue-associated micro- 
biota can be traits of Escherichia and Enterococcus. The 
virulence-associated functions such as pore formation by 
hemolysin E [7], bacterial attachment to mammalian 
cells via AIDA-I adhesin-like protein [8], cytotoxicity to 
eukaryotic cells by cytolethal toxins [9] and cytotoxic 
necrotizing factor 1 [10], increased motility (outer mem- 
brane usher protein, fimbrial-like protein, type 1 fimbrial 
protein, fimbrial chaperone protein) [11], and invasion 
of epithelial cells via invasion B [12] can be attributed to 
Enterobacteriaceae, particularly, Escherichia. The proteo- 
lytic cleavage by coccolysin also known as gelatinase can 
be attributed to Enterococcus [13,14]. Importantly, gelati- 
nase GelE of Enterococcus faecalis has been shown to 
degrade collagen, a critical component for anastomotic 
tissue healing [15,16]. We further performed enrichment 
analysis for KEGG metabolic pathways to determine 
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Figure 3 Comparative analysis of bacterial abundance in anastomotic tissues at PODO and POD6. The analysis is presented at phylum, 
order, class, family, and genus levels. 



what metabolic changes occur at post-anastomotic tis- 
sues versus stool. The analysis indicates there are sub- 
stantially more metabolic changes for tissue-associated 
microbiota than for the microbiota in the luminal con- 
tents, as expected based on changes in microbial com- 
munity structure (Figure 4). Furthermore, the analysis of 
post-anastomotic tissues predicts virulence-associated 
metabolic pathways such as lipopolysaccharide biosyn- 
thesis, bacterial secretion system, and shigellosis, which 
are not represented in the stool microbiota. The meta- 
bolic pathways were linked to the OTUs responsible for 
the functions (Figure 5). The analysis demonstrated a 
shift in the primary contributors to the metabolic func- 
tions such as genus Lactobacillus, order Clostridiales, 
and genus Ruminococcus for pre-anastomotic tissue ver- 
sus family Enterobacteriaceae, genus Enterococcus, and 
genus Bacteroides for post-anastomotic tissue (Figure 5). 

Discussion 

Advances in non-culture-based analyses of the intestinal 
microbiome have demonstrated that environmental in- 
fluences (e.g., diet, antibiotic therapy, age, colonization 
history, geography) have profound effects on microbiota 
community structure, membership, and function [18], 
which are region specific within the gastrointestinal tract 
[4,5]. Although intestinal surgery in the form of resec- 
tion and anastomosis has been practiced for decades, 
until now, there has been no systematic assessment of 
the bacterial diversity and functional potential associated 
with this surgery. It has been observed in germ-free and 
antibiotic-treated animals that disruption of the normal 
microbiota results in a weaker intestinal anastomosis as 
measured by burst strength (i.e., the pressure at which 



an anastomosis will burst) [19,20]. In these studies, they 
also observed a loss of collagen, a major determinant 
of anastomotic healing and burst strength. Conversely, 
empiric antibiotic treatment in other animal models pre- 
vents anastomotic leak when pathogenic bacteria or 
pathologic conditions are present [21]. This latter obser- 
vation suggests that tissue injury alone, and the physio- 
logic stress it incurs, may significantly influence the 
microbiome. 

Our results indicate that intestinal resection and anas- 
tomosis induce a change in the microbiota associated 
with intestinal tissues, which is not observed in the 
luminal contents. Our findings also predict that tissue- 
associated microbial communities express an adhesive 
phenotype, which may explain, in part, their tropism to- 
ward injured tissues. Previous work from our laboratory 
has indicated that soluble products released during sur- 
gical injury can act as both chemoattractants and 'cues' 
that induce a phenotype shift among pathogenic bacteria 
such as P. aeruginosa [22,23]. In the present study, anas- 
tomotic tissues may select for microbes that express 
enhanced virulence (i.e., E. faecalis) and thus may have 
the potential to complicate anastomotic healing. A crit- 
ical loss of the cytoprotective microbiota (i.e., Firmicutes 
and Bacteroidetes) on anastomotic tissues may further 
facilitate colonization and invasion by disease-associated 
microbiota, again with the potential to complicate heal- 
ing. The precise explanation by which the process of 
surgery in the present study shifted the intestinal micro- 
biota remains to be clarified. The surgery performed in 
the present study included a large abdominal incision, 
colon resection, anastomosis, and wound closure, all of 
which has been shown to result in a catabolic stress and 
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Table 1 Predicted upregulated functions in tissue-associated microbiota after anastomosis 


KEGG 
orthology 


Description 


Pathway/function 


Primary 
contributor 


p value 


K03276 


UDP-glucose/galactose:(glucosyl)LPS alpha-1, 
2-glucosyl/galactosyltransferase [EG2.4.1 .-] 


Lipopolysaccharide biosynthesis 


Enterobacteriaceae 


0.02 


K1 1 139 


Hemolysin E 


Pore-forming toxin [7] 


Enterobacteriaceae 


0.02 


K03765 


Transcriptional activator of cad operon 


The cad operon encodes a system for 
neutralization of low extracellular pH 


Enterobacteriaceae 


0.02 


K12678 


AIDA-I adhesin-like protein 


Autotransporter family porin, mediates bacterial 
attachment to mammalian cells [8] 


Enterobacteriaceae 


0.02 


K09925 


Hypothetical protein 






0.02 


K11014 


Cytolethal distending toxin subunit B (CdtB) 


Initiates a eukaryotic cell cycle block at the G2 
stage prior to mitosis, CdtB potentiates a 

II I ■ II 1 1 1 1 rr\~\ 

cascade leading to cell cycle block [9] 


Enterobacteriaceae 


0.02 


K1 1015 


Cytolethal distending toxin subunit C (CdtC) 


CdtC and CdtA function as dimeric subunits, 
which bind CdtB and delivers it to the mammalian 
cell interior 


Enterobacteriaceae 


0.02 


K11013 


Cytolethal distending toxin subunit A (CdtA) 


Same as above 


Enterobacteriaceae 


0.02 


K1 1 264 


Methylmalonyl-CoA decarboxylase 

rcr - -A 1 1/111 
LbL.4. 1 . 1 .4 1 J 


Propanoate metabolism 


Enterobacteriaceae 


0.02 


K1 3244 


c-di-GMP-specific phosphodiesterase 
LtL.3. 1 .4.jzJ 


Hydrolases acting on ester bonds 


Enterobacteriaceae 


0.02 


K03459 


Formate transporter (focB) 


Electrochemical potential-driven transporters 


Enterobacteriaceae 


0.02 


K07354 


Outer membrane usher protein (sfmD) 


Bacterial motility proteins [11] 


Enterobacteriaceae 


0.02 


K07355 


Fimbrial-like protein (sfmF) 


Bacterial motility proteins 


Enterobacteriaceae 


0.02 


K07352 


Type 1 fimbrial protein (sfmA) 


Bacterial motility proteins 


Enterobacteriaceae 


0.02 


K07353 


Fimbrial chaperone protein 


Bacterial motility proteins 


Enterobacteriaceae 


0.02 


K12985 


(Galactosyl)LPS 1 ,2-glucosyltransferase 
[EC2.4.1-] 


Lipopolysaccharide biosynthesis proteins 


Enterobacteriaceae 


0.02 


K10708 


Fructoselysine 6-phosphate deglycase 
[EC3.5.-,] 


Hydrolases acting on carbon-nitrogen bonds, other 
than peptide bonds 


Enterobacteriaceae 


0.02 


K11008 


Cytotoxic necrotizing factor 1 (cnfl) 


Blocks cell cycle G2/M transition in uroepithelial cells [10] 


Enterobacteriaceae 


0.02 


K1 1 920 


AraC family transcriptional regulator (envY) 


DNA-binding transcriptional regulator 


Enterobacteriaceae 


0.02 


K13285 


Invasin B 


Bacterial invasion of epithelial cells [12,17] 


Enterobacteriaceae 


0.03 


K08605 


Coccolysin [EC34.24.30] (gelE) 


Metallopeptidases [13,14] 


Enterococcus 


0.03 



immune activation as a compensatory mechanism to 
deal with the process of healing the wounds [24,25]. 
There likely exist multiple mechanisms that drive the 
compositional and functional changes in the microbiota 
including redistribution of nutrients away from the gut 
to healing tissues, immune activation, and changes in lu- 
minal and mucosal oxygen status [26]. It is also possible 
that the microbiota 'sense' a dramatic change in host health 
status and undergo compositional and functional changes 
as a mechanisms of adaptation until recovery to homeosta- 
sis is complete. Further work is needed to determine pre- 
cisely which of the factors predominate in this model. 

It is noteworthy that most of the current antibiotic 
regimens used as prophylaxis agents for intestinal sur- 
gery are based on culture-based analyses of expelled 
stool. The limitation of this approach is that antibiotics 
are not based on a tissue-specific microbiome. As such, 



we do not know which tissue microbiota are associated 
with healing versus non-healing. This may be critically 
important to prevent complications such as anastomotic 
leaks, which still persist today at unacceptable rates even 
under the most expert of care [27]. While antibiotic use 
has been demonstrated to be highly effective in prevent- 
ing anastomotic leak, the type, duration, and microbial 
targets of currendy recommended antibiotics have remained 
empirical and invalidated. Data from the present study 
would suggest that direct genomic and functional analyses 
of the tissue-associated microbiota at anastomotic sites may 
provide a more directed approach to antibiotic use in anas- 
tomotic surgery by understanding their effect on the micro- 
bial community structure, membership, and function. 

The observed 500-fold increase in Enterococcus in post- 
operative anastomotic tissues is interesting. Although en- 
terococci are commensal inhabitants of the mammal gut, 
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Figure 4 (See legend on next page.) 
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(See figure on previous page.) 

Figure 4 Enriched metabolic pathways for pre- and post-anastomic tissue and stool. Pathway enrichment for KEGG metabolic pathways 
using HUMAnN followed by statistical comparative analysis using LEfSe was performed to determine differential enrichment between PODO and 
P0D6 for tissue and stool. Red color indicates that the KEGG metabolic pathway is enriched in comparison to the corresponding POD for either 
tissue or stool. Tissue and stool were not directly compared with each other, so common metabolic pathways do not indicate equivalent levels 
of enrichment for both. 



they have been observed to be regularly associated with 
human infections [28] . In terms of intestinal anastomoses, 
it has been reported that E. faecalis is highly prevalent at 
leaking anastomotic tissues, although its role as a causa- 
tive agent in leakage remains unknown [29,30]. Ente- 
rococcus has a high adherence affinity to extracellular 
cellular matrix proteins such as fibronectin, laminin, and 
various different types of collagens [31-33] including colla- 
gen IV [34]. What role, if any, Enterococcus plays in the 
etiopathogenesis of anastomotic leak will require further 
study. However, the fact that its abundance is increased to 
such a high degree raises the possibility that it may com- 
plicate healing given that it is known to confer a high degree 
of pathogenicity when it predominates. In fact, PICRUSt 
analysis [6] in the current work predicted the capacity of 
anastomotic tissues-associated microflora to possess cocco- 
lysin known as a gelatinase GelE of £ faecalis that is able to 
degrade collagen. Similarly, Escherichia genus also displayed 
a dramatic increase in relative abundance (i.e., 200-fold) at 
anastomotic tissues which is another commensal capable of 
confer a pathogenic phenotype [35]. Further work will be 
necessary to independentiy assess the virulence potential of 
strains isolated from anastomotic tissues versus luminal con- 
tents in both cases. 

Given that the normal microbiota provide pathogen 
colonization resistance [36-38], there may be value in 
understanding how the microbiota assemble from the 
lumen to anastomotic tissues over the course of anasto- 
motic healing. For example, our results demonstrated 
that there is a postoperative decrease in Ruminococca- 
ceae, a family of Firmicutes that belong to obligate 
anaerobes. It has been recently shown that the families 
Ruminococcaceae and Lachnospiraceae (Clostridium clus- 
ter XIV group) colonize mucosal folds in the proximal 
murine colon [5]. Taken together, our data indicate that a 
significant membership of the normal microbiota may be 
lost at anastomotic tissues, the consequence of which is 
unknown. The family Lactobacillaceae was also reduced 
at anastomotic tissues and its known cytoprotective action 
[39-41] to accelerate healing may play a role in leak when 
they are critically decreased and additional physiologic 
stress occurs. 

Although taxonomic profiling by the 16S rRNA ana- 
lysis can be very informative, it does not provide any 
functional context to the bacterial communities. Metage- 
nomics and metatranscriptomics would provide greater 



detail on their potential and actual functions yet are still 
technically challenging in terms of analysis. However, 
using PICRUSt [6], we observed predicted global changes 
in metabolic and virulence functions of bacterial com- 
munities among microbiota associated with anasto- 
motic tissues. These findings provide a path forward in 
which we can apply this technique to models of anas- 
tomotic leak and determine its predictive capacity and 
its ability to inform a microbial mechanism for anasto- 
motic leakage. 

Conclusion 

An anastomotic injury appears to significantly alter the 
tissue-associated microbiota without affecting the lu- 
minal microbiota. A more detailed analysis of the tissues 
associated microbiota over the full course of anastomotic 
healing could lead to novel strategies for a leak preven- 
tion that may include allowing for a balance of microbial 
community structure, membership, and function. 

Methods 

Rat model of anastomosis 

The Institute for Animal Care and Use Committee at 
the University of Chicago approved all animal experi- 
ments. Adult, male Wistar rats 250-300 g (Charles River 
Laboratory, Chicago, IL, USA) were used for all experi- 
ments and were allowed unrestricted access to rat chow 
and tap water. Prior to surgery, rats were sedated with 
40-80 mg/kg ketamine, 5-10 mg/kg xylazine via an in- 
traperitoneal injection. Using aseptic technique, rats 
were subjected to an abdominal incision (laparotomy) 
and a 0.5-cm distal colon resection followed by an end- 
to-end rectosigmoid reconnection (anastomosis) using 
12 interrupted 6-0 non-absorbable sutures. Colon lu- 
minal content was collected for DNA isolation. The 
resected piece of colon was washed with saline and then 
swabbed for tissue associated microbiota. In all cases, 
anastomotic integrity was confirmed immediately after 
construction by injecting saline into the rectum and 
observing that the anastomotic was water tight. On post- 
operative day 6 all rats were sacrificed and the luminal 
contents adjacent to anastomosis were collected. For 
tissue-associated microbiota, anastomotic tissues were 
removed, and mechanically homogenized, and then DNA 
was extracted. 
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Figure 5 Contribution of taxonomic groups to the enriched metabolic pathways for pre- and post-anastomotic tissue. There is a shift in 
the primary contributors to the metabolic functions from taxa such as genus Lactobacillus, order Clostridials, and genus Ruminococcus for pre-anastomotic 
tissue to taxa such as family Enterobacteriaceae, genus Enterococcus, and genus Bacteroides for post-anastomotic tissue. Values are scaled to a 0-1 range for 
each KEGG metabolic pathway (ko) individually, with 1 indicating the highest contributing OTU and 0 indicating no contribution by the OTU. 



16S rRNA analysis of bacterial composition 

The following are the four groups of samples in our 
studies: (1) PODO tissue for which the resected piece of 
colon was washed with saline and then swabbed for 
tissue associated microbiota; (2) PODO luminal bacteria 
representing luminal content from a resected piece of 
colon; (3) POD6 tissue for which tissues at the site of 
anastomosis were collected, washed with saline, and then 



swabbed; (4) POD6 luminal bacteria representing luminal 
content from collected piece of anastomotic site of the 
colon. Twenty rats were included in the experiment. 
Samples with sequences containing low confidence base 
calls (phred < 20) were eliminated. Samples were analyzed 
by 16S rRNA V4 iTAG amplicon sequencing analysis. 
Following joining of paired 2 x 150 bp MiSeq read 
sequences (Illumina, San Diego, CA, USA) containing 
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low confidence base calls (phred < 20) were dropped, 
and the remaining sequences clustered into 97% simi- 
larity operational taxonomic units (OTUs) using uclust 
[42,43]. After dropping clusters composed of a single 
sequence, OTUs were assigned taxonomic annotations by 
aligning representative members of each cluster to the 
greengenes 12 04 database using the RDP classifier 
[44,45]. Sequencing depth among samples was normal- 
ized by rarefaction to a depth of 3,000 sequences per 
sample. Along with clustering, taxonomic assignment, 
and rarefaction, PCoA plots were derived from weighted 
unifrac distance matrices using the QIIME vl.6.0 software 
suite [46]. Significant differences in the abundance of 
taxonomic groups between treatments were calculated 
with Student's t test followed by the Benjamini-Hochberg 
correction for multiple comparisons using an alpha of 
0.05 [47]. Taxa tree was generated using GraPhlAn. 

Functional analysis of bacterial composition 

We predicted the functional composition and abundance 
of the microbial community using PICRUSt on the 
OTUs derived from the 16S rRNA analysis [6]. Differen- 
tial functional abundance between postoperative day 0 
and day 6 for stool contents and anastomosis tissue was 
performed using the Wilcoxon rank sum statistical test, 
and p values were Bonferroni corrected for multiple test- 
ing. The metagenomic predictions produced by PICRUSt 
were further analyzed by HUMAnN [48] to provide 
KEGG metabolic pathway coverage, and LEfSe [49] 
was used for comparative analysis of the differentially 
enriched KEGG metabolic pathways between postoper- 
ative day 0 and day 6 for stool contents and anasto- 
mosis tissue samples. OTU contribution for enriched 
metabolic pathways was calculated by summing the abun- 
dance contribution of each OTU for each KEGG orthology 
comprising a metabolic pathway, and then normalizing by 
the number of samples. As there are hundreds of OTUs, a 
threshold was used to select just the top contributing 
OTUs for inclusion in Figure 5. The values were scaled to 
a 0-1 range for each individual metabolic pathway, with 1 
indicating the highest contributing OTU and 0 indicating 
no contribution by that OTU. 
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POD: postoperative day; KEGG: Kyoto Encyclopedia of Genes and Genomes (a 
database resource for understanding high-level functions and utilities of the 
biological system); PICRUSt: Phylogenetic Investigation of Communities by 
Reconstruction of Unobserved States (a bioinformatics software package 
designed to predict metagenome functional content from marker gene (e.g., 
16S rRNA) surveys and full genomes); OTUs: operational taxonomic units. 
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